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ABSTRACT 

The continuous temporal coverage and high photometric precision afforded by space observatories 
has opened up new opportunities for the study of variabihty processes in young steUar cluster members. 
Of particular interest is the phenomenon of deuterium-burning pulsation in brown dwarfs and very- 
low-mass stars, whose existence on 1-4 hours timescales has been proposed but not yet borne out 
by observations. To investigate short-timescale variability in young, low-mass objects, we carried out 
high-precision, high-cadence time series monitoring with the Warm Spitzer mission on 14 low mass 
stars and brown dwarfs in the ^-^3 Myr a Orionis cluster. The flux in many of our raw light curves is 
strongly correlated with sub-pixel position and can vary systematically as much as 10%. We present 
a new approach to disentangle true stellar variability from this "pixel- phase effect," which is more 
pronounced in Warm Spitzer observations as compared to the cryogenic mission. The light curves 
after correction reveal that most of the sample is devoid of variability down to the few-millimagnitude 
level, on the minute to day timescales probed. However, one exceptional brown dwarf displays erratic 
brightness changes at the 10-15% level, suggestive of variable obscuration by dusty material. The 
uninterrupted 24- hour datastream and sub-1% photometric precision enables limits on pulsation in 
the near-infrared. If this phenomenon is present in our light curves, then its amplitude must lie below 
2-3 millimagnitudes. In addition, we present three field eclipsing binaries and one pulsator for which 
optical ground-based data is also available. 

Subject headings: open clusters and associations: individual {a Orionis) — stars: low-mass, brown 
dwarfs — stars: variables: general — techniques: photometric 
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1. INTRODUCTION 

Photometric monitoring of brown dwarfs (BDs) and 
very low mass stars (VLMSs) in young clusters provides 
insight into the dynamic processes affecting such objects 
at a few million years of age, including accretion, mag- 
netic effects, and star-disk interaction. It has long been 
known that T Tauri stars and their higher mass counter- 
parts exhibit optical brightness fluctuations of typically 
a few tenths of a magnitude, but in many cases more 
than a magnitude on timescales of hours to days. How- 
ever, few observations have explored the short-timescale 
regime (i.e., seconds to hours) in the lowest mass ob- 
jects. This is in part due to the difficulty of obtain- 
ing from ground-based facilities suitably long time se- 
ries data with high cadence and minimal interruption. 
Nevertheless young BDs and VLMSs may display sig- 
nificant unexplored variability. Objects with active ac- 
cretion can exhibit luminosity fiuctuations related to in- 
falling gas on second to minute timescales, owing to the 
instability of the shock position with respect to the stella r 
photosphere (|Koldoba et al.l I2008L lOrlando et al.ll2010( ). 
BDs and VLMSs whose surfaces are not obscured by in- 
falling material are expected to display a different type of 
instability- pulsation fueled by central deuterium burn- 
ing (e.g., ;Palia fc Baraffc 200l hereafter PB05). The pe- 
riods for this phenomenon are predicted to be in the 1-4 
hour range for objects from ~0.02 to O.IMq. Empirical 
verification of these instability theories through detec- 
tion of short-timescale aperiodic or periodic variability 
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presents a new opportunity to probe the properties of 
young low-mass objects. 

A few previous claims of short-period variability 
in the optical have reported flux changes at t h e few 
to five percent level (e.g., iScholz &: Eisloffell 12004 
Bailer-Jones fc Mundtl 1200 IL iZapatero Osorio et all 



20031) . However, our own attempts to identify pulsation 



from th e ground in a number of very-low-mass a Ori 
objects (|Codv fc HiUen brandllMol hereafter CHIO) re- 
sulted in more stringent limits on pulsation and shorter 
timescale periodicities in the optical: If present, /-band 
amplitudes must be less than ~0.01 magnitudes in the 
brown dwarf sample, and an order of magnitude lower 
among the very low mass stars. Thus further efforts in 
the search for and characterization of short-timescale 
variability require even more sensitive observations. 

Space telescopes offer a chance for deeper variability 
searches since the lack of atmosphere minimizes system- 
atic errors in photometry, affording signal-to-noise ratios 
close to the poisson limit. They also fulfill the need for 
dense and continuous time sampling by staring at a sin- 
gle patch of sky for extended periods of time without the 
inconveniences of weather, daytime interruption, or syn- 
optic scheduling. Additional progress may be made by 
observing in the infrared. While this band is not a tradi- 
tionally favored band for photometric time series work, it 
has several advantages for the detection of low-amplitude 
variability in BDs. Because of their cool temperatures, 
BDs are brightest at wavelengths just longward of 1 /xm 
and thus should be amenable to relatively high signal-to- 
noise photometry in the near to mid-infrared. Our opti- 
cal observations revealed that approximately 85% of low- 
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mass cluster members display variability at the 1-10% 
level, which we attributed to primarily rotational modu- 
lation of spots and variable accretion. The amplitude of 
brightness fluctuations produced by these mechanisms is 
expected to decrease with wavelength (e.g. iFrasca et al.l 
[2009i) , thereby reducing confusion between pulsation and 
other sources of variability. Thus while the amplitude 
range and wavelength dependence of pulsation are un- 
known (the linear stability theory of PB05 predicts only 
periods, as a function of mass), the lower temperature 
contrast between any magnetic spots or accretion flows 
and the photosphere may enhance the detection proba- 
bility in the infrared. 

We have used the Spitzer Space Tele scope warm 
mission Infrared Array Camera (IRAC; iFazid 12004 
IWcr ncr 2003) to observe a sample of 14 brown dwarfs 
and low-mass stars in the ^3 Myr a Orionis clus- 
ter, several of which had been claimed previously as 
short-period variables ( e.g.. iBailer-Jones fc Mundt 2001, 
iZapatero Osorio et al.l 120031 iSchoirfcEisloffel ,2004) . 
The reported timescales (2-5 hours) are too short to be 
explained by rotational modulation, suggesting that at 
least some of these objects are good accretion instabil- 
ity or pulsation candidates. With data at the 3.6 and 
4.5 /Ltm wavelengths, the pixel-phase effect (i.e., oscil- 
lations in the measured flux introduced by uncorrected 
intrapixel sensitivity variations) is more pronounced now 
than in the cryogenic mission. In this paper, we discuss 
an alternative method relative to the commonly adopted 
approach to remove it. We present light curves sampling 
timescales from one minute to one day and discuss the 
prospects for pulsation. Finally, we identify several field 
eclipsing binary and pulsating variables for which opti- 
cal light curves are also available from our ground-based 
study, and compare the behavior at both wavelengths. 

2. TARGET SELECTION AND PROPERTIES 

Many very-low-mass young cluster members are now 
catalogued and thus available for time series monitor- 
ing. Since Spitzer/IRAC has two fields that are each 
5.22' across, only a few clusters in the 1-10 Myr range 
contain enough known very low mass members to en- 
able monitoring of more than one or two BDs simulta- 
neously. Among these are IC 348 and a Orionis. We 
chose the latter for the present study with Spitzer and 
the former for investigation with the Hubble Space Tele- 
scope (results to be presented in a future paper) . Within 
a Ori, the objects S Ori 3 1, S Ori 45 have been c l aimed 
as short-term variables (iBailer- Jones fc Mundtl I2001L 
IZapatero Osorio et al.l I2003D . providing two promising 
targets for our pulsation search. Since the region around 
them is encompassed by our previously published optical 
fields (CHIO), we used the same list of candidate clus- 
ter members to select targets. We experimented with 
the position of the 3.6 fim field as well as its orientation 
with respect to the 4.5 /xm field, whose center is offset by 
^^6.7', to optimize the pointings. 

In addition, we aimed to select targets with a high 
probability of exhibiting pulsation, based on luminosity 
and temperature consistent with PB05's predicted po- 
sition of the pulsation instability strip. To assess H-R 
diagram positions, we assembled a set of 50 previously 
confirmed a Ori members with available spectral types, 
primarily from .Barrado v Navascucs et al.. (,2003. ). Tem- 



peratures were estimated using the intermediat e grav- 
ity temperature scale derived by iLuhman et al.l ([2003) , 
which accounts for the lower gravity of young objects 
compared to field dwarfs and is appropriate for the young 
objects studied here. In addition, t hey hav e been cali- 
brated for consistency with the Bar affe et al.l fl998.) low- 
mass evolutionary models, on which the pulsation in- 
stability strip from PB05 is based. For a few objects 
without prior spectral types, we obtained low-resolution 
(i? ^1400) spectra from the Double Spectrograph on the 
200-inch Hale Telescope at Palomar Observatory. Spec- 
tral types were estimated by comparison with data taken 
with the same set-up for ^3- Myr low-mass IC 348 m em- 
bers previously classified by ILuhman et al.l (|2003D . as 
well as ~1 Myr Tau rus and ^5 M yr Upper Scorpius 
members observed bv lSlesnick et al.l |2006a.b). Our new 
spectra, along with those of a number of other objects 
in the a Ori cluster will be presented in a forthcoming 
paper. We adopt uncertainties of 100 K, equivalent to 
just under one spectral subclass. 

Luminosities of a Ori members are dependent upon the 
estimated distance to the cluster. This value has often 
been taken to be 350lgQ" pc, based on the Hipparcos par- 
allax of cr Ori AB itself. However, iSherrv et aH ()2008f ) 
showed that a distance of 440^30 is more consistent with 
main sequence fitti ng to observations of cluster A stars. 
iJeffries et a l. (2006) pointed out that what has tradition- 
ally been considered the a Ori cluster is in fact likely a 
superposition of two kinematically distinct groups with 
different radial velocities, ages, and distances. They pro- 
pose that one of the populations corresponds to the Orion 
OBla and OBlb association subgroups, while the other 
is associated with the star a Ori its elf. With the s e con - 
siderations in mind, we adopt the ISherrv et all (|2008D 
distance but for completeness we also explore (in §5.1) 
the effect of the smaller value on our computed lumi- 
nosities and positions on the H-R diagram. The result- 
ing distance moduli, m-AI, are 8.21±0.15 and 7.72 ±0.65 
magnitudes. Extinction tow ard a Ori is relatively low, 
and we adopted Aj=0.044 (|Barrado v Navascues et all 

Hoop). 

Final luminosities were determined with J-band 
magnitudes fro m Barrado v Navascues et al.l ([2003), 
iCaballerol (2008), and Bcia r et all ((19991) ■ Both the J 
and / bands arc generally favored for their relative lack 
of contamination from accretion and disk excess. How- 
ever, bolometric corrections in J have the additional ad- 
vantage of b eing less sensit ive to color and surface grav- 
ity age (e.g.,|Luhman"1999'). We adopted the bolometric 
corrections used in Kraus & Hillenbrand (2007), and also 
used those of lCaballero et al.l ()2007[ ) to check that the re- 
sults were relatively insensitive to the form of the correc- 
tions as a function of color and spectral type; we adopt 
their value of 0.15 magnitudes as a typical uncertainty. 

We show the computed locations of our late-type ob- 
jects on the H-R diagram with respect to the theoreti- 
cal pulsation instability strip in Fig. 1, for both possi- 
ble distance modulus values. Uncertainties in luminos- 
ity include photometric and bolometric correction errors. 
However, the true errors are dominated by the system- 
atic uncertainty in the distance to the cluster, as shown 
in the figure. 

Additional systematics may be introduced by the 
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Fig. 1. — The instability strip for pulsating brown dwarfs and 
very low mass stars is shown in blue, along with isochrones 
from [Baraffe et al.l l|2003l ) and a sample of spectroscopically con- 
firmed low-mass members of t he a Ori cluster drawn from 
IBarrado v Navascues et al.l II2003I V In the top panel we assumed 
a distance of 440 pc, while in the bottom we changed this to the 
alternate value of 350 pc. Targets observed with Spitzcr IRAC ap- 
pear in red. A number of these objects are on or near the predicted 
instability strip, suggesting that they might exhibit pulsation. 

choice of band used to calculate the luminosity. We per- 
formed a comparison test of luminosities derived from 
the / band for those objects in our Spitzer sample with 
available /-band photometry. There is an approximately 
uniform discrepancy of ^0.35 dex between luminosities 
derived from the J-band magnitudes, versus the /-band 
magnitudes. One might conclude that the J-band mag- 
nitudes include contributions from circumstellar disks, 
but in fact the J-band luminosities are fainter. Such a 
discrepancy may be caused by the unknown difference 
between the dwarf-like bolometric corrections adopted 
here and those that account for the lower surface gravi- 
ties of young objects. We retain the luminosities as de- 
rived from the J band but in §5.1 compare results from 
both optical bands. 

The final Spitzer sample- shown in red in Fig. 1- in- 
cluded five confirmed and two candidate BDs in a Orio- 
nis, with three in the 3.6 fim field and four in the 4.5^m 
field. In addition, we observed serendipitously seven 
other known a Ori cluster members in the fields which 



likely are too massive (e.g., >0.1 Mq) to exhibit pulsa- 
tion but are nonetheless valuable targets for investigation 
of other types of young star variability. This brought the 
total sample to 14 objects- six in the 3.6 fiui field, and 
eight in the 4.5 fim field. Fewer objects were placed in 
the 3.6 /im field because of scheduling constraints on the 
required orientation. Details on each target are provided 
in Table 1 (see Table 1 in CHIO for more details, includ- 
ing coordinates as well as 2MASS identifications). All 
except S Ori 31 and S Ori 53 are spectroscopically con- 
firmed members of the a Ori cluster; both have colors 
and spectral type consistent with low-mass a Ori mem- 
bership, while the form er also has a proper motion consis- 
tent with membership (jLodieu et al.ll2009l) . In addition, 
we consider object SWW 188 a new spectroscopically 
confirmed member since it exhibits weak Na I absorp- 
tion indicative of low surface gravity in the low- resolution 
spectra that we obtained. 

3. WARM SPITZER OBSERVATIONS AND DATA 
REDUCTION 

Prior predictions and limits on the amplitudes (<0.01 
magnitudes) and timescales (~l-4 hours) of pulsation 
guided our observational setup. The ability to detect 
light curve periodicities at a particular amplitude (A) 
depends on both the photometric noise level (cr) as well 
as the total number of data points (A''). When identi- 
fying a periodic signal in a periodogram (e.g., §5), the 
signal-to-noise ratio in frequency space is roughly equal 
to A^/N /(2a) and must be larger than --4.0 for 99.9% 
certainty (see the discussion in §5.1 of CHIO). We set 
a target of several millimagnitudes for the minimum de- 
tectable periodic amplitude in all objects apart from the 
faintest two BDs. In addition, data had to be taken fre- 
quently enough to probe periodicities on timescales close 
to an hour. Accordingly, observations were carried out 
over a 24-hour period from 22 to 23 October 2009 (Astro- 
nomical Observation Request key 35146240 and program 
identification 60169). Exposure times were 23.6 seconds 
each, resulting in a cadence of ^30 seconds, and a total 
of 2730 data points. 

Observations at both wavelengths take place simulta- 
neously, with one of the fields in each of the 3.6 /Ltm 
and 4.5 /im cameras. Therefore, we collected data only 
in a single band for each of our targets. The orienta- 
tion of the two fields is shown in Fig. 2, and the centers 
were R.A.=05'^38"23.3", decl.=-02°40'29" (3.6 /im) and 
R.A.=05'^38'"26.4^ decl.=-02°47'13" (4.5 /xm). The po- 
sition angle was -94.7° east of north for both fields. Since 
our aim was to produce photometric time series with as 
high a precision as possible, we elected not to dither. 
Keeping the positions of all sources fixed within a single 
pixel reduces the effect of flux variations introduced by 
pixel-to-pixel sensitivity differences not fully corrected by 
flatfielding, although intrapixel sensitivity variation (the 
"pixel-phase effect") remains an issue and is addressed 
below. 

For data acquired from the Spitzer /ISlkC camera, all 
basic calibrations are performed via pipeline, as ex- 
plained in the handboolfl. As of version 18.12.0, the 
IRAC pipeline provided images at several different stages 
of processing, from raw unreduced frames to final phot- 
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TABLE 1 
Basic target data 



Object 


/ 


3.6 


4.5 


SpT 


Ref 


Optical variability 


4771-41 


12.95±0.02 


- 


8.84±0.02 


K5 


5 


Aperiodic: RMS=0.23 mags 


SWW40 


14.18±0.03 


- 


11.61±0.01 


M3 


5 


Periodic: 4.47d, 0.013 mags 


S Ori J053817.8-024050 


15.00±0.04 


11.68±0.01 




M4 


5 


Periodic: 2.41d, 0.008 mags 


SWW188 


15.06±0.03 




12.61±0.01 


M2 


5 




S Ori J053823. 6-024132 


15.13±0.04 


12.17±0.05 




M4 


5 


Periodic: 1.71d, 0.017 mags 


S Ori J053833. 9-024508 


16.15±0.04 




12.52±0.03 


M4 


5 


Aperiodic: RMS=0.06 mags 


S Ori J053826. 8-022846 


16.17±0.04 


12.71±0.03 




M5 


5 




S Ori J053825. 4-024241 


16.96±0.04 


12.96±0.03 




M6 


2 


Aperiodic: RMS=0.16 mags 


S Ori J053826. 1-024041 


17.05±0.04 


13.65±0.01 




M6 


2 




S Ori J053829. 0-024847 


17.06±0.05 




12.91±0.03 


M6 


3 




S Ori 27 


17.22±0.05 


13.13±0.01 




M7 


1 




S Ori 31 


17.46±0.04 




13.67±0.02 


M7 


1 




S Ori 45 


20.03±0.09 




15.05±0.05 


M8.5 


1 


Periodic: 0.3d, 0.034 mags 


S Ori 53 


20.31±0.09 




17.5±0.4 


M9 


4 





Note. — Wc list the 14 confirmed and candidate a Orionis cluster members observed with Spitzer^ in order of optical brightness, /-band 
magnitudes are taken from CHIO. 3.6 and 4.5 /im band photometry is the median value determined over our light curves, with conservative 
uncertainties including systematic errors due to poor knowledge of intra-pixel sensitivity distributions as well as intrinsic variability. Values listed 
in the optical variability column are either the RMS spread of aperiodic light curves over a ~2-week period, or the period and amplitude of periodic 

light curves. 

References- (1) IBarrado v Navascues et al.l II2003I) : (2) ICaballero et al.l II2006I) : (3) ICaballero et all 120071) : (4) IBarrado v Navascues et al.l II2001I) : 
(5) this work. 

ometry ready data. However, at the time of writing there 
were stih a number of problems resulting from the tran- 
sition to the Warm Spitzer mission. Standard bias and 
dark subtraction, flatfielding, linearity and flux calibra- 
tions have been applied to create the basic calibrated 
data (BCD) files. Further corrections including auto- 
mated removal of cosmic rays and the column pull-down 
effect have been performed to create a set of corrected 
BCDs (CBCDs). Since these procedures were fine-tuned 
to cryogenic mission data, they left numerous column 
pull-down artifacts as well as a residual bias pattern in 
our data. Therefore, we elected to carry out the final set 
of reductions manually, starting with the BCDs. 

Because there are no laboratory-generated bias frames 
corresponding to warm mission conditions, we retained 
the bias subtraction applied by the pipeline and mod- 
eled the remaining uncorrected pattern. Fortunately the 
residuals largely consist of vertical bands in which bright- 
ness remained relatively constant throughout our obser- 
vations. A procedure to median stack all images for each 
channel, mask out the objects, and reset each column to 
its mean value was performed by S. Carey (2010, private 
communication). Subtraction of the resulting vertical 
striped bias correction image from all BCDs effectively 
removed the residual patterns. 

The column pull-down effect, in which counts are re- 
duced throughout columns with bright (>35,000 DN) 
sources, was also not fully corrected for in the pipeline. 
Unlike in cryogenic mission data, flux values associated 
with pull-down now differ above and below the source, 
in addition to following an approximately exponential 
trend as a function of y position on the array. We were 
provided an updated pulldown correction code (D. Pal- 
adini 2010, private communication), which satisfactorily 
modeled and removed this effect. 



4. PHOTOMETRY ROUTINE 




Fig. 2. — The Spitzer IRAC 3.6 and 4.5 fim fields are overlaid on 
a Palomar Observatory Sky Survey 2 (POSS2) red image obtained 
from the Digitized Sky Survey (DSS)^. A portion of our 20x20' 
/-band ground-based campaign with the CTIO 1.0-meter telescope 
as published by OHIO is also shown for reference (dashed region). 

We performed aperture photometry on our 14 target 
objects using a variety of aperture sizes and sky annulus 
widths and radii for background subtraction. The IRAC 
camera has a somewhat undersampled point spread func- 
tion (PSF), with a full- width at haff maximum (FWHM) 
size of ~1.4 pixels. As a result, most of the flux from 
each object is concentrated in the central pixel. Inaccu- 
rate aperture centering can thus lead to erroneous bright- 



4.1. Aperture photometry 



^ http:/ /archive. stsci.edu/cgi-bin/dss_form 



ness fluctuations in the resulting light curve. We deter- 
mined moment centroid positions by calculating position- 
weighted flux averages within a four pixel radius. Points 
for which the centroid algorithm failed due to a cosmic 
ray or other bad pixel effect were omitted from the data 
(this comprised ^-^3% or less of the photometry). Aper- 
tures with radii from two to four pixels were placed at 
these centroids, sky annuli of various sizes from 2 to 9 
pixels were used, and the enclosed sky-subtracted flux 
was determined with the IRAF phot task. We adopted 
the aperture resulting in the lowest RMS light curve 
spread, which was 2 pixels for most targets. We con- 
verted the data to magnitudes by incorporating the pub- 
lished IRAC zero point values, aperture corrections, and 
location-dependent array response provided by the hand- 
book. Average 3.6 and 4.5 fim band magnitudes are 
listed in Table 1. 

Even with careful placement of apertures, many of the 
light curves contained deviations beyond the expected 
white noise level that were not characteristic of the un- 
derlying stellar variability. Points with particularly large 
flux suggested cosmic ray hits within the stellar PSF. 
These occurrences appear random and uncorrelated, and 
thus unlikely to represent real short-term astrophysical 
behavior. Since we did not dither, it was not possible to 
remove these without binning images or datapoints. We 
elected instead to filter erroneous flux values directly out 
of the light curves with a 3-sigma clipping algorithm. 

4.2. Pixel phase correction 

Our first pass at the photometry also revealed that 
most objects suffer from the well known IRAC pixel- 
phase effect: although target positions were restricted 
to a single pixel, movement of the centroid within the 
pixel introduces position-correlated flux changes of up 
to 1 0% due to response variations within individual pix - 
els (|Morales-Calder6n et ahl 120061 iDeming et al.ll2011[ ). 
The X and y centroid positions executed not only several 
small jumps, but also an oscillatory motion with period 
~60 minutes due to the subtle effect of a thermal cycling 
battery heater on Spitzer pointing. As a result, most of 
the light curves from channel 1 exhibited periodic fluc- 
tuations of up to ^4% amplitude, along with additional 
systematics of up to 10%, or 0.1 magnitudes. We display 
a typical example of x and y trends as a function of time 
in Fig. 3. The effect is about half as large in channel 2 
but still signiflcant enough to require removal in many of 
the light curves. 

The Warm Spitzer mission guide presents a method 
to correct these effects by providing a model of the sen- 
sitivity variations within a pixel. However, the model 
was derived from observations of a single bright star and 
does not account well for differences in the response pat- 
terns of different pixels. We found that the proposed 
correction algorithm was not adequate for removing the 
pixel-phase related noise from our light curves. Typi- 
cal signal-to-noise ratios were 55-60% of that estimated 
based on the Poisson limit, whereas previous work with 
warm Spitzer data suggests that we should be a ble to 
achieve upwards of 75-80% (jPeming et al.l [20Tll) . On 
the other hand, subtraction of a median-fit trend from 
each light curve confirmed that the white noise level did 
indeed reach a level consistent with these predictions. 

To recover the additional ~20% in S/N, we explored 
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Fig. 3. — X and Y pixel centroid positions of one of our targets 
(S Ori J053817.8-024050) as a function of time. Since telescope 
pointing affects all targets in a similar fashion, both the short- 
timescale (~0.04d) oscillations and the more systematic trends are 
typical of the centroid behavior of other observed objects as well. 
Outlier points indicate where the centroiding algorithm has failed 
(e.g., because of a cosmic ray hit or other artifact). 

several methods for removing noise due to the pixel-phase 
effect. The failure of the model based on a single bright 
star implied that the spatial response differs significantly 
from pixel to pixel. Therefore, we attempted to fit each 
object's flux with polynomials as a function of x and y 
position. Unfortunately this approach proved problem- 
atic for several reasons. First, the pointing during our 
Astronomical Observation Request (AOR) traced out a 
region in x-y space that was neither homogeneous nor 
large compared to the pixel size (e.g.. Fig. 4). Rather, 
small pointing jumps led to centroid positions occupying 
three somewhat discrete areas of the pixel. In addition, 
we were concerned that intrinsic variability of our young 
cluster sources could complicate the fitting process. 

Plots of flux versus x, y, and phase (distance from a 
flxed point near the center of the pixel) did not exhibit 
tight trends, suggesting that accurate removal of system- 
atic effects would not be feasible. As a result, we opted 
to flt a gaussian functional form to each object's spa- 
tial flux distribution. The Warm Spitzer guidqj suggests 
that a double gaussian function (i.e., sum of gaussians in 
the X and y directions) is the best-fitting pixel sensitiv- 
ity model based on bright star data. However, because 
of the incomplete spatial coverage within each pixel, we 
suspected that a single gaussian with adjustable center 
would work as well. Our adopted pixel sensitivity model 
thus consisted of four free parameters: 

/;s^F(.-a^-^af+iv-yof}/2<7^ + Fq, 

^ http:/ /ssc. spitzer.edu/irac/warmfeatures 
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Fig. 4.— Variations in the flux of object S Ori J053817.8-024050 
measured within a single pixel. The pixel is centered at [0,0] and 
extends to ±0.5 in the x and y directions. Only a portion of the 
pixel is depicted here. 

where AF is the height of the gaussian function, x and y 
are the positions with respect to the center of the pixel, 
Xo and yo are the offsets of the peak flux response from 
the center of the pixel, and a is the width of the gaussian. 
Fq is a constant determined so that the function averages 
to 1.0 over the entire pixel. 

To identify the best-fitting pixel-phase function we cre- 
ated a script to iterate through a reasonable range in 
the four parameters, perform the pixel phase correction 
based on each set, and assess the presence of pixel-phase 
noise in the resulting light curve. This assessment was 
performed by generating a periodogram in the range of 
frequency space where the pixel-phase oscillation domi- 
nates: 21.5-25 cycles/day (corresponding to periods of 
hours, and unfortunately very close to the pulsa- 
tion signature that we seek) . It is here that a large peak 
is seen in the periodograms of raw light curves (Fig. 5). 
We present as the "corrected" light curve the one for 
which the integrated periodogram in this region is mini- 
mized. In some cases (S Ori 31, S Ori J053833.9-024508, 
SWW 188, S Ori 45, and S Ori 53), the initial peri- 
odogram did not display a peak associated with the pixel- 
phase effect, and so we did not apply any correction. 
Since the correction process only targets a small region 
of frequency space in the periodogram, it should preserve 
variability that is intrinsic to the objects, if present. 

We emphasize that we have chosen the symmetric 
gaussian pixel-phase model out of convenience and lack 
of knowledge of the underlying distribution; the true 
pixel sensitivity function is likely to be much more com- 



plicated fe.g.. iBallard et ani2010D . The presented light 
curves may thus have systematic inaccuracies. In ad- 
dition, since the correction process removes only varia- 
tion on the known ~ 1-hour period of the thermal oscilla- 
tion, it is not obvious as to whether variation on longer 
timescales is intrinsic to the sources or undercorrection of 
the pixel-phase effect. We caution that any Warm Spitzer 
studies attempting to assess variability whose precise na- 
ture (i.e., light curve shape) is not known in advance will 
face this issue. 

5. VARIABILITY SEARCH RESULTS 

Our main aim in analyzing the light curves of low- 
mass a Ori cluster members is to search for periodici- 
ties on the 1-4-hour timescales predicted for deuterium- 
burning inst ability. We pro duced fourier transform pe- 
riodograms (lD eemin3ll975l ) for both the raw and pixel- 
phase-corrected light curves; all are presented in Fig. 5. 
Since the observations were continuous over a 24-hour 
period at 30-second cadence, we are sensitive to periodic 
variability on timescales from one minute to one day. In 
addition, the periodogram does not suffer from aliasing, 
so true signals are relatively easy to identify if they rise 
high enough above the noise baseline. In most cases, the 
periodograms display a relatively uniform mean from fre- 
quencies at a few cycles per day (cd~^) out to the Nyquist 
limit at 1440 cd~^. This white noise level depends on the 
magnitude of the source and ranged from 0.001 to 0.004 
magnitudes in the periodogram. 

Examination of the periodograms revealed that the 
pixel-phase correction process substantially lowered the 
noise level, enabling better sensitivity to periodicities 
outside the 1-1.1-hour range of the pixel-phase oscilla- 
tion. The two exceptions were SOri 27 and 4771-41. The 
former was centered near the edge of two pixels, making a 
fit to the spatial distributions difficult without resorting 
to a more complex non-gaussian function. Object 4771- 
41 is exceedingly bright, and residual variability seen in 
the final light curve may be a figment of the correction 
process. 

The majority of periodograms are relatively feature- 
less, refiecting minimal variability in the input light 
curves. In the low- frequency region from one to sev- 
eral cd~^, many of the periodograms steadily rise in a 
"1//" fashio n indicative of sy stematic or "red" noise 
trends (e.g., iPress et al.l Il992j ) close to the maximum 
timescale of one day. However, apart from one object 
(SOri J053825. 4-024241), no potential periodic signals 
stand out high enough above the baseline. Here we have 
used a criterion of 4-cr, equivalent to 99.9% certainty, 
as explained in CHIO. Since SOri J053825.4-024241 dis- 
plays intrinsic aperiodic variability at the 10% level, and 
the putative 3.5 cd^^ (or ^9.6h) signal barely exceeds the 
detection threshold at S/N~4, it is unclear as to whether 
this is a true periodicity. 

5.1. Prospects for pulsation 

The lack of periodic signals in the 1-4-hour range sug- 
gests that none of the a Ori cluster members in our 
sample exhibits deuterium-powered pulsation at a level 
above several millimagnitudes. However, the strength 
of this conclusion depends on the likelihood that one 
or more targets fall on the PB05's predicted pulsation 
instability strip. If we assume a distance of 350 pc 
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Fig. 5. — Light curves and periodograms for all targets, in order of decreasing optical brightness. Object identifications are listed above 
each light curve, and the Spitzer band (3.6 or 4.5 /im) is noted in the y-axis label. Black indicates the raw light curve and periodogram, 
whereas red shows them after correction for intra-pixel sensitivity effects. Objects with no red points did not require correction. Periodogram 
frequencies are given in cycles per day. Insets show the same periodograms zoomed in to the low-frequency range where the signature of 
the pixel- phase oscillation is visible (~22-24 cd~^). 
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for the a Ori cluster, then all seven BDs in the sam- 
ple (S Ori J053825.4-024241, S Ori J053826.1-024041, 
S Ori 31, S Ori J053829.0-024847, S Ori 53, S Ori 27) 
may be on the instability strip, to within the uncertain- 
ties. If we instead adopt a distance of 440 pc, then the 
VLMS S Ori J053826.8-022846 becomes an additional 
candidate, whereas the position of S Ori 45 falls slightly 
off the strip. Thus one would naively expect that a hand- 
ful of our targets have temperatures and luminosities 
consistent with those required for pulsational instabil- 
ity. Nevertheless, the significant size of the measurement 
uncertainties compared with the width of the strip must 
be taken into account. 

To estimate the chance that in fact none of our sam- 
ple have H-R diagram positions overlapping the instabil- 
ity strip, we adopt temperature-luminosity probability 
distributions for each object. We take the distributions 
to be two-dimensional asymmetric gaussians, normalized 
and centered at the adopted luminosities and tempera- 
tures. The gaussian widths are given by the associated 
1-(T uncertainties. The position of each target then corre- 
sponds to a probability that it is susceptible to pulsation, 
which we determine by integrating its distribution over 
the entire region of the instability strip. For objects on or 
very close to the strip, this value is ~20-25%, whereas for 
the higher mass stars far from the strip it is close to zero. 
The probability that the position of a given object does 
not overlap with the instability strip then ranges from 
75% to 100%. The product of these values over all tar- 
gets provides an estimate of the chance that no pulsators 
would be present in our sample. For luminosities derived 
from J-band magnitudes and a cluster distance of 440 pc, 
this probability is 19%- small but certainly not negligi- 
ble. Since systematic errors can significantly affect the 
result, we have also recalculated this value using both I- 
band-derived luminosities and the alternate distance of 
350 pc. The different combinations yield probabilities 
from 23% to 32%. Turning these numbers around, there 
is a ~70-80% chance that at least one object should ex- 
hibit pulsation based on its position within the instability 
strip, assuming that the theoretical calculations under- 
pinning it (PB05) do not suffer from gross systematic 
errors. 

In addition, the expectation value for the number of 
objects lying directly on the strip lies between one and 
two, depending on the choice of distance and J or / band 
magnitudes. Therefore if pulsation is operating at an ob- 
servable level, we are likely to detect at least one instance 
of it. The fact that we also did not detect short-timescale 
variability in our larger ground-based sample (CHIO) of 
BDs and VLMSs suggests that the lack of observable pul- 
sation is a repeatable result and not the product of sam- 
ple selection effects. Nevertheless, we caution that the 
small number of low-mass objects observed here leaves 
open the possibility that no true pulsation candidates 
were among our sample. 

Assuming this is not the case, statistically we expect 
at least one object in our sample to be susceptible to 
pulsation. If so, the amplitude of this phenomenon ap- 
pears to be too low to be observed in the infrared. To 
quantify the detection limit, we have fit power laws to 
each periodogram, tracing out the maximum amplitude 
level as a function of frequency. These curves, of form 
A/ [f + B) -\- C for frequency / and constants A, B, and 



C, mark the highest amplitude signal observed for each 
object, whether real or noise. We display them in Fig. 
6, where the procedure has been carried out on both the 
raw and pixel-phased corrected data. 

For the majority of objects, we have detected no pe- 
riodicities in the pulsation frequency range with am- 
plitudes greater than several millimagnitudes. Brown 
dwarfs S Ori 45 and S Ori 53 are exceptions, as they have 
higher limits (0.005 and 0.04 magnitudes respectively in 
the 4.5 fim band) owing to their faintness and corre- 
spondingly high noise levels in both the light curves and 
periodograms. In addition, brown dwarf S Ori J053825.4- 
024241 has a higher limit for pulsation (0.004-0.007 mag- 
nitudes in the 3.6 /zm band, depending on frequency) 
since it displays substantial intrinsic variability. The rest 
of our targets have maximum amplitudes in the peri- 
odogram of at most 0.002 to 0.003 magnitudes. This 
represents the threshold above which we detect no pe- 
riodicities. We conclude that if deuterium-burning pul- 
sation is present in any of our sources (apart from the 
three exceptions noted above), then its amplitude must 
be below this level. 

5.2. SOriJ053825. 4-02^241: a high- amplitude variable 
brown dwarf 

Among our sample, the substellar a Ori member 
SOriJ053825.4-024241 stands out as the lone target 
highly variable on timescales less than 24 hours. With a 
3.6 fim band RMS of 0.035 magnitudes, this object has a 
peak-to-peak amplitude of 0.15 magnitudes. It displays 
variations about four times as large in the /-band, based 
on our longer timescal e ground-bas ed dataset (CHIO). 
Other studies ( Caballcr o et al.|[2006[ ) have indicated that 
SOriJ053825. 4-024241 is actively accreting and has a disk 
iHcrnandcz et al. 2007). 

No previous infrared studies of brown dwarfs have un- 
covered aperiodic variability on such short timescales. 
However, variability of young stars at Spitzer wave- 
lengths or of brown dwarfs in general with these am- 
plitudes and on longer time scales is not unprecedented. 
The Yo ung Stellar Object Variability (YSOVAR) project 
(jMorales-Caldero n 2011) campaign on young Orion Neb- 
ula Cluster stars (masses >0.1 Mq) has also found sub- 
stantial erratic variability in the 3.6 and 4.5 /zm bands. 
Assessment of their data has shown that the aperi- 
odic variables among the sample known to harbor disks 
display a range of variability RMS values centered on 
^0.03 magnitudes in the 3.6 /im band (Morales-Calderon 
2011, private communication). Similar amplitude distri- 
butions were obtained using existing multi-epoch data 
with limited c a dence in T aurus and Chamaeleo n I by 
ILuhman et all (|2008D and ILuhman et all (|2010D . The 
typical RMS of a few hundredths of a magnitude is 
quite consistent with the valu e that we have measured 
for S Ori J053825.4-024241. iMorales-Calde^ (|20TT[ ) 
discuss the possible causes of the mid-infrared variabil- 
ity and surmise that many of their variables may be 
explained by variable obscuration by overdense regions 
in the inner disk, while others are caused by intrinsic 
changes in the inner disk emission itself. Either of these 
scenarios may apply to SOriJ053825. 4-024241. In any 
case, hot accretion gas is likely too faint at infrared wave- 
lengths to serve as the source of variability for this object. 

To further explore the behavior of this BD on different 
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Fig. 6. — Limits on pulsation detection in the periodograms based on raw light curves. Solid red curves show the limits for objects on 
or near the instability strip, whereas grey dashed curves are for objects not expected to exhibit pulsation. The curve for S Ori 53 is off 
the chart at a uniform amplitude of 0.04 magnitudes. The dotted lines indicate the region of frequency space where we expect pulsation to 
occur (i.e., 1-4-hour periods). The left panel displays data before pixel-phase correction, and the right shows the final data after removal 
of pixel-phase oscillation trends. For the few objects that did not require these corrections, the curves derived directly from raw data are 
repeated. 



timescales, we have performed an autocorrelation analy- 
sis. In addition to displaying quasi-periodicity patterns 
not picked up by the periodogram, it is useful in assess- 
ing the timescale on which the variability mechanism re- 
mains coherent. We have calculated an autocorrelation 
function based on the S Ori J053825. 4-024241 hght curve 
using both a standard, "biased" formula, as well as one 
that corrects for the finite data length. The standard 
autocorrelation function (ACF) is given by: 

N-t/At 

^w^iT^ E mw+t/At), 

where L{j) are the light curve points, Ai is the time 
spacing between datapoints (which must be uniform) , N 
is the total number of points, and the A(Q) factor in 
front in included so that at a time lag of zero, the ACF 
is completely correlated (^(0) = 1). 

To account for the fact that fewer points are avail- 
able to calculate the ACF at longer lag times {t > 0.5), 
we have produced another version- the "unbiased" ACF- 
in which the this roughly linear effect (~ iV — t/At)) 
has been divided out. In both cases, we have computed 
the autocorrelation via fourier transform of the power 
spectrum (as specified by the Wiener-Khinchin theorem; 
IWieneri[T930llKhinchinl 11931 . since this is both faster 
and less prone to numerical inaccuracies. 

Both versions of the ACF are plotted in Fig. 7. We find 
that the light curve is well correlated up to timescales of 
~0.15d, or 3.6h. At longer timescales, it also shows sig- 
nificant correlation due to the overall trend seen in the 
light curve; this is illustrated by the two peaks at ~0.43d 
and ~0.9d (the latter primarily in the unbiased ACF). 
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Fig. 7.— Autocorrelation functions for SOriJ053825.4-024241. 
The solid line is the standard ACF, whereas the dashed line is the 
version that has been corrected for finite data length. 

We conclude that the variability mechanism is physically 
coherent on timescales of at least a few hours. The hy- 
pothesis of variable obscuration in association with the 
disk is qualitatively consistent if the scale of dumpiness 
and location of dust is such that fluctuations would pass 
by the face of the BD in several hours. 

6. COMPARISON WITH OPTICAL DATA 

Previously we observed a region including both Spitzer 
flelds, using the CTIO 1.0-meter telescope (CHIO). High 
cadence (every 10 minutes) data in the Cousins / band 
and lower cadence (twice per night) data in R was ac- 
quired over runs of 12 and 11 consecutive nights, respec- 
tively in 2007 and 2008. We identified a number of vari- 
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able objects, both in the a Ori cluster and background 
field. Although the overall time baseline of the Spitzer 
observations is short compared to that of the ground- 
based campaign, we have searched for common variabil- 
ity in the two datasets. 

A total of seven variable a Ori cluster members from 
the ground-based campaign fall in the fields of the Spitzer 
observations, as noted in Table 1. In the 4.5 /im field, 
S Ori J053833.9-024508 and 4771-41 were identified as 
aperiodic variables in the ground-based photometry. In 
addition, the BD S Ori 45 was identified as being periodic 
in the /-band, with a period of 7.2 hours and amplitude 
0.034 magnitudes, whereas VLMS SWW40 was found to 
have a period of 4.47 days and amplitude 0.013 mag- 
nitudes. In the 3.6 fim field, the BD S Ori J053825.4- 
024241 was identified as an aperiodic variable, as dis- 
cussed in §5.2. Two additional variables were found 
be periodic in the / band: S Ori J053823.6 (P=1.7d; 
^=0.017 mags) and S Ori J053817.8-024050 (P=2.4d; 
A=0.008 mags). 

For those ground-based variables with brightness fluc- 
tuations on timescales longer than a day, we do not neces- 
sarily expect to observe variability in our shorter Spitzer 
dataset. Indeed, we do not recover periodic variability 
at great than the 1% level in any of the ground-based 
periodic variables. In addition to the shorter time base- 
line, it is possible that the non-simultaneity of observa- 
tions and the different wavelengths make rotational spot 
modulation- the primary explanation for periodic vari- 
ability in young VLMSs and BDs- unobservable in our 
light curves. 

Several of the previously identified aperiodic variables, 
on the other hand, do appear to be variable at infrared 
wavelengths. The BD S Ori J053825.4-024241 displays 
relatively high amplitude erratic fiuctuations (see §5.2). 
Object 4771-41 shows residual variability after correc- 
tion for the pixel-phase effect (light curve RMS of ^0.01 
versus ~0.001 magnitudes), and S Ori J053833.9-024508 
also displays variability at a significantly higher level 
than predicted by signal-to-noise estimates (light curve 
RMS of -0.05 versus <0.01 magnitudes). The RMS val- 
ues in the Spitzer bands are similar to those found in the 
optical for S Ori J053825.4-024241 and S Ori J053833.9- 
024508, whereas they are roughly an order of magnitude 
lower for 4771-41. Thus the light curve of this latter ob- 
ject may exhibit residual pixel-phase effects, as opposed 
to real variability. However, for the other two aperiodic 
variables, the rough correspondence of RMS amplitudes 
in both the optical and infrared suggests that the vari- 
ability mechanism may be relatively insensitive to wave- 
length. 

Interestingly, object S Ori J053829.0-024847 displays 
substantial variability at 4.5 /im (an 0.06-magnitude drift 
over 24 hours), whereas it did not appear variable in 
our ground-based dataset. We suspect that the variabil- 
ity mechanism in this case was dormant during the op- 
tical observations, although its photometry could have 
been affected by a nearby neighbor on the array. Since 
this object exhibits an infra red excess (iHernand ez et al.l 
[2007L ICaballero et al.ll2007l) . there is an additional pos- 
sibility that the variability is associated with the disk 
and thus only visible in the near-infrared and at longer 
wavelengths. 

In addition to the recovery of aperiodic variability in 



our a Ori cluster sample, we also re-identify a number 
of eclipsing binaries; further details on these field objects 
are provided in the appendix. 

7. CONCLUSIONS 

We have presented high-cadence infrared light curves 
of 14 low-mass a Orionis cluster members based on Warm 
Spitzer observations. The excellent precision of our pho- 
tometry led to limits of 0.002-0.003 magnitudes on the 
amplitude of any deuterium-burning pulsation in these 
objects in the 3.6 and 4.5 fim bands. This result is con- 
sistent with our ground-based /-band findings, which re- 
vealed no periodicities with amplitudes greater than 0.01 
magnitudes and timescales shorter than 7 hours among 
low-mass a Ori cluster members. In this work we have 
reduced the amplitude limit by an order of magnitude, 
albeit on a smaller sample of BDs and VLMSs. No- 
tably, we also find little other variability on the < 24- 
hour timescales probed by the data. The main exception 
was brown dwarf S Ori J053825. 4-024241, which displays 
brightness variations of up to 0.1 magnitudes over the 
course of a day. The similarity of amplitudes in the in- 
frared and optical suggests that obscuration by dust ma- 
terial in the surrounding disk provides a better explana- 
tion than does variable light scattering or accretion (e.g., 
Ivan Boekel et ahllMl . We propose that a general lack 
of variability among young, low-mass cluster members 
may in fact be useful for future studies at infrared wave- 
lengths, such as searches for planets around young BDs 
and VLMSs. Both transit detection and radial veloc- 
ity measurements benefit from low levels of spot-related 
activity on short timescales. 

We also emphasize that production of light curves 
devoid of the pixel-phase and other detector effects is 
difficult at present with Warm Spitzer data. While 
previous work has successfully identified low-amplitude 
planetary transit si gnatures in Spitzer light curves (e.g., 
iDeming et al.|[20TlT ). the transit event represents only a 
small portion of these time series and hence systematic 
trends can be taken into account. When the entire light 
curve is instead the subject of interest, and the form of 
variability is unknown in advance, systematics are more 
difficult to model and remove. For future high-precision 
photometric time series work, we recommend further ex- 
ploration of the sensitivity distribution within individ- 
ual pixels, perhaps through even higher cadences that 
might provide more datapoints over a given time and 
thus greater spatial coverage within individual pixels. 

This work is based on observations made with the 
Spitzer Space Telescope, which is operated by the Jet 
Propulsion Laboratory, California Institute of Technol- 
ogy under a contract with NASA. We acknowledge sup- 
port from NASA under contract 1382589 administered 
through JPL / Caltech. We thank the referee, Kevin Luh- 
man, for helpful comments. We also express appreci- 
ation to Sean Carey and Roberta Paladini for help in 
correcting pull-down and bias residuals, as well as Maria 
Morales-Calderon for extensive advice on Warm Spitzer 
reductions. 
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APPENDIX 

In addition to examining the light curves of the 14 

(T Ori cluster member targets, we also searched the en- 
tire 3.6 and 4.5 /im fields for serendipitous foreground 
and background variables. After producing light curves 
for all point sources with magnitudes less than ^ 19.0, 
we assessed their RMS spread as a function of bright- 
ness. Objects lying more than three standard deviations 
above the median trend were flagged as possible vari- 
ables. We visually examined their light curves and disre- 
garded those whose brightness fluctuations were clearly 
caused by pixel sensitivity effects. Four objects (other 
than BD 053825.4-024241; §5.2) displayed conspicuous 
variability by these criteria; their light curves are pre- 
sented in Fig. 8. For consistency with the other presented 
light curves, we show both the time series and their pcri- 
odograms. We list the estimated period, which often does 
not correspond to the largest periodogram peak since this 
analysis method is relatively insensitive to the presence 
of secondary eclipses. 

All four stars were also identified as variables in our I- 
band ground-based datasct (CHIO); therefore, we refer to 
them by the same nomenclature. We have not rigorously 
fit eclipse profiles or other models to the data but present 
estimates (~ 10-20% accuracy) of their main parameters 
here: 

CTIO J05381870-0246582 is an eclipsing binary 



system with an /-band depth of ^0.45 magnitudes, and 
4.5 /im depth of at least 1.2 magnitudes. The most likely 
period is ~11.8 hours, or 5.9 hours if all of the eclipses 
are primary (the data are too noisy to to distinguish dif- 
ferent depths in subsequent eclipses). 

CTIO J05382129-0240318 also appears to be an 
eclipsing binary, with period ^9.6 hours. This period 
is fully consistent with our ground-based data, for which 
we unfortunately reported an erroneous value (4.6 days 
instead of 9.5 hours). The 3.6 /im depth (> 1.3 mag- 
nitudes) is significantly deeper than the 7-band depth 
(~0.35 magnitudes). 

2MASS J05382188-0241039 exhibits a slightly 
asymmetric periodic profile, reminiscent of an RR Lyrae 
star. Its period of 11. 8h is also consistent with this type 
of pulsator. Since the timescale is so close to half a day, 
aliasing caused us to misidentify and report a l.Od pe- 
riod for the ground-based data. The 3.6 /xm peak-to-peak 
amplitude is ^0.25 magnitudes, whereas the value at / 
band is just over 0.6 magnitudes. 

2MASS J05381949-0241224 also displays the char- 
acteristic shape of a close eclipsing binary, although there 
is slight decrease in its peak amplitude over 24 hours 
which may be attributed to systematic pixel sensitiv- 
ity effects. The period is 2.8 or 5.6 hours, depending 
on whether alternating brightness dips are secondary 
eclipses. At ~0.5 magnitudes, the peak-to-peak ampli- 
tude at 3.6 /um is about 20% smaller than that in the / 
band. 
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Fig. 8. — Field variable stars. Light curves (top) and periodograms (bottom) are as in Fig. 5; estimated periods are marked near 
corresponding frequency peaks. 
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